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ABSTRACT 

We have obtained 17 moderate-resolution (~ 2.5 A) optical spectra of the Galactic 
X-ray Nova Velorum 1993 in quiescence with the Keck-II telescope. Cross-correlation 
with the spectra of late-type dwarfs (especially K7-M0) yields excellent radial 
velocities for the secondary star. The orbital period (P) is 0.285206 ± 0.0000014 d, 
and the semiamplitude (K2) is 475.4 ± 5.9 km s~ 1 . Our derived mass function, 
/(Mi) = PK\j2'KG = 3.17 ± 0.12 Mq, is close to the conventional absolute limiting 
mass for a neutron star (~ 3.0-3.2 Mq) — but if the orbital inclination i <^ 80° (given 
the absences of eclipses), then Mi ;> 4.2 — 4.4 Mq for nominal secondary-star masses 
of 0.5 Mq (M0) to 0.65 M (K6). Even if the secondary is quite undermassive (e.g., 
M2 = 0.3 Mq), we derive Mi J> 3.9 Mq. The primary star is therefore almost certainly 
a black hole rather than a neutron star. 

Fits to the wings of the double-peaked Ha emission line yield approximate radial 
velocities for the compact primary. The velocity curve has a semiamplitude (K\) of 
65.3 ± 7.0 km s , but with a phase offset by 237° (rather than 180°) from that of 
the secondary star. Under the assumption that the observed semiamplitude reflects 
the motion of the primary (despite this offset, which is common), the mass ratio 
q = M 2 /Mi = Ki/K 2 = 0.137 ± 0.015, and hence for M 2 = 0.5-0.65 M we derive 
Mi = 3.64-4.74 Mq. The constraints from q and the mass function yield Mi = 4.4 Mq 
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and i w 78° if we use a normal K7-K8 secondary (M2 ~ 0.6 Mq). Indeed, consistency 
cannot be achieved for M2 <^ 0.59 Mq if the maximum inclination estimate (80°) is 
correct. An adopted mass Mi ~ 4.4 Mq is significantly below the typical value of 
~ 7 Mq found for black holes in other low-mass X-ray binaries. 

Keck observations of MXB 1659-29 (V2134 Oph) in quiescence reveal a probable 
optical counterpart at R = 23.6 ± 0.4 mag. This object recently went into a new X-ray 
and optical (R pa 18 mag) outburst, after about 21 years of inactivity. 

Subject headings: binaries: close — black hole physics — X-rays: stars — stars: 
individual (Nova Vel 1993, MXB 1659-29) — novae, cataclysmic variables 

1. INTRODUCTION 

An important class of binary systems has been identified in which a low-mass secondary 
(companion) star orbits a probable black hole; see White, Nagase, & Parmar (1995), van Paradijs 
& McClintock (1995), and Tanaka & Lewin (1995) for extensive reviews. In all cases they were 
first observed in outburst as "X-ray novae" (or "soft X-ray transients"). Their X-ray spectra are 
generally characterized by a prominent, "soft" thermal component (kT ~ 1 keV) as well as a 
"hard" power-law tail extending to very high energies (0.1-1 MeV). During outburst, the radiation 
from these and other low- mass X-ray binaries (LMXBs; the "low mass" refers to the secondary 
star) is emitted predominantly by the accretion disk surrounding the primary star. 

A lower limit to the mass of the primary (Mi) in a given X-ray transient can be measured 
when it returns to quiescence. At that time, light from the secondary star contributes significantly 
to (or even dominates) the visible spectrum, and the secondary's radial-velocity curve can be 
determined with a series of time-resolved spectra (see Cowley 1992 for a review). The orbital 
period (P) and semiamplitude (K2) of the secondary yield the mass function of the primary, 
/(Mi) = PK^jlixG = Mf sin 3 i/(Mi + M2) 2 , where i is the inclination of the orbital plane to our 
line of sight. Clearly, /(Mi) provides an absolute lower limit to the mass of the primary; only if 
M2 = and i = 90° is Mi = /(Mi). It is generally acknowledged that if the primary is dark and 
has /(Mi) ^ 3.2 M Q , it is probably a black hole, since the theoretical upper limit to the mass of a 
normal neutron star is ~ 3.0-3.2 Mq (Friedman, Ipser, & Parker 1986; but see Friedman & Ipser 
1987, as well as Bahcall, Lynn, & Selipsky 1990, for conditions that might allow neutron stars to 
exceed this nominal limit). 

The six best examples of LMXBs whose mass function exceeds (or is close to) ~ 3 Mq , along 
with the derived values of /(Mi) (in units of Mq) and references, are as follows in chronological 
order of the mass-function measurement: A0620-00 = V616 Mon (3.18 ± 0.16, McClintock & 
Remillard 1986; 2.72 ± 0.06, Marsh, Robinson, k Wood 1994; 2.91 ± 0.08, Orosz et al. 1994), 
GS 1124-68 = Nova Mus 1991 (3.1 ± 0.4, Remillard, McClintock, & Bailyn 1992; 3.01 ± 0.15, 
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Orosz et al. 1996), GS 2023+338 = V404 Cyg (6.26 ± 0.31, Casares, Charles, k Naylor 1992; 
6.08 ± 0.06, Casares & Charles 1994), GRO J1655-40 = Nova Sco 1994 (3.16 ± 0.15, Bailyn et 
al. 1995; 3.24 ± 0.09, Orosz & Bailyn 1997), GS 2000+25 = Nova Vul 1988 (5.02 ± 0.47, Casares, 
Charles, & Marsh 1995; 4.97 + 0.10, Filippenko, Matheson, & Barth 1995a, slightly revised to 
5.01 ± 0.15 [not ±0.12] by Harlaftis, Home, & Filippenko 1996), and Nova Oph 1977 (4.0 ± 0.8, 
Remillard et al. 1996; 4.86 ± 0.13, Filippenko et al. 1997, slightly revised to 4.65 ± 0.21 by 
Harlaftis et al. 1997). Here we add a seventh object to this list: Nova Vel 1993 = GRS 1009-45, 
with a derived mass function of 3.17 ± 0.12 Mq. 

Nova Vel 1993 was discovered on 12 September 1993 with the WATCH all-sky monitor 
aboard Granat (Lapshov, Sazanov, & Sunyaev 1993; Lapshov et al. 1994) and with BATSE on the 
Compton Gamma-Ray Observatory (Harmon et al. 1993). Its spectrum exhibited an "ultrasoft" 
hump at low energies (<; 1 keV) and a power-law tail out to at least 100 keV (Kaniovsky, Borozdin, 
&: Sunyaev 1993), typical of X-ray binaries in which the compact object is a black hole. Two 
months later (17 November), Delia Valle & Benetti (1993) discovered a blue optical counterpart 
at V ~ 14.6 mag, but reasonable estimates suggest that the magnitude at the time of outburst 
was V = 13.8 ± 0.3 (Delia Valle et al. 1997). Optical photometry conducted by Bailyn & Orosz 
(1995) about half a year after the primary outburst showed the presence of a secondary outburst 
and several mini-outbursts, again reminiscent of black-hole X-ray novae. Delia Valle et al. (1997) 
suggested an orbital period of about 4 hours, but a more reliable period of 6.86 ± 0.12 hours was 
obtained by Shahbaz et al. (1996). The spectral type of the secondary star in the binary system 
was estimated to be late-G/early-K by Shahbaz et al. (1996), and later than G5-K0 by Delia 
Valle et al. (1997). 

On 1998 January 25 (UT dates are used throughout this paper), we obtained several -R-band 
images of the field of Nova Vel 1993 with the Low Resolution Imaging Spectrometer (LRIS; Oke et 
al. 1995) at the Cassegrain focus of the Keck-II telescope. As can be seen in Figure 1, which shows 
a subset of one image (seeing 0.65"), the nova was in quiescence by this time, with R = 21.2 ± 0.2 
mag.^ To determine whether Nova Vel 1993 should be considered a good dynamical black-hole 
candidate, we decided to obtain a radial- velocity curve with LRIS. Our group had already obtained 
excellent results in this way for the black-hole candidates GRO J0422+32 (Filippenko, Matheson, 
& Ho 1995b; Harlaftis et al. 1999), GS 2000+25 (Filippenko et al. 1995a; Harlaftis et al. 1996), 
and Nova Oph 1977 (Filippenko et al. 1997; Harlaftis et al. 1997), all of which are comparably 
faint. 



2 We adopt the magnitudes of comparison stars quoted in Table 1 of Delia Valle et al. (1997). There appears to be 
a numbering mismatch, or errors in the photometry, of some stars in Table 1 of Shahbaz et al. (1996): for example, 
their Stars 2, 5, and 7 should have comparable magnitudes, yet they are listed as being very different. 
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2. OBSERVATIONS AND REDUCTIONS 

Nova Vel 1993 was observed with LRIS in 1998 during the nights of January 25, February 
1, March 5-6, and May 2, as well as in 1999 during the night of January 21. A journal of useful 
observations is given in Table 1. (The spectra obtained on 1998 May 2, and a few spectra on other 
nights, were of marginal quality and are not considered here.) Given the object's far southerly 
declination (—45°), and the restrictive southwest azimuth limit on Keck-II (<^ 185°), observing 
was restricted to ~ 1.5 hours per night; hence, a number of different observing runs separated by 
a range of intervals was needed to avoid serious aliasing of the orbital period. Conditions were 
always clear, and the seeing was about 1.2", reasonably good considering the high airmass (~ 2.4). 
Typical exposure times were 1100 s. 

The long slit of width 1" was oriented at a position angle (PA) of 160° for all observations 
except that of 1998 January 25. This was close to the parallactic angle at the time of observation, 
thereby reducing differences in the relative amount of light lost at different wavelengths (Filippenko 
1982). At such a PA, the slit went directly through Star A, which is much brighter than the 
quiescent nova and only ~ 1.6" SE of it. (Delia Valle et al. 1997 incorrectly state that Star A is 
SW of the nova.) However, since the direction of atmospheric dispersion (i.e., the parallactic angle) 
at the time of observation coincided with the nova/Star-A orientation, the degree of contamination 
from Star A was nearly independent of wavelength. The slit also partially intersected Stars B 
and C (Fig. 1), but their light did not contaminate that of the nova; see Figure 2, which plots 
the intensity of light along the slit at the wavelength of Ha (thereby accentuating the nova's 
contribution). 

We used a Tektronix 2048 x 2048 pixel CCD with a scale of 0.215" pixel" 1 (0.43" per binned 
pixel in the spatial direction). The 1200 grooves mm -1 grating, blazed at 7500 A, resulted in 
a wavelength range of ~ 5650-6950 A, and a full- width at half- maximum (FWHM) spectral 
resolution of ~ 2.5 A (~ 120 km s _1 ) with the 1" slit, essentially identical to what we used in our 
previous studies of black-hole X-ray novae. Thirteen velocity standards with spectral types in 
the range G5 V-M2 V were observed with the same setup, given the estimated classification of 
the secondary star in Nova Vel 1993 (Shahbaz et al. 1996; Delia Valle et al. 1997). Also, spectra 
of some sdF stars (Oke & Gunn 1983) were obtained for flux calibration and removal of telluric 
absorption lines. 

Cosmic rays were eliminated from the two-dimensional spectra through comparison of pairs 
of consecutive exposures. The two-dimensional spectra were bias-subtracted and flattened in the 
usual manner. The wavelength scale was determined from polynomial fits to the positions of 
emission lines in spectra of Hg-Ne lamps obtained with the telescope at (or near) the position 
of each object. To ensure accurate wavelength calibration, final corrections (~ 0.0-0.3 A) to the 
wavelength solution were obtained from night-sky emission lines in the spectra of the nova. 
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We used the APALL task in IRAF^ to optimally extract (Home 1986) one-dimensional, 
sky-subtracted spectra of Nova Vel 1993. However, to minimize contamination from Star A, we 
extracted only a few rows of the CCD centered on the position of the nova; as shown in Figure 2, 
these were rows 3-5 from the center of Star A, in the left wing of its spatial profile, corresponding 
to a displacement of 1.3-2.1". To further remove contamination, we also extracted the same rows 
(3-5) in the right wing of Star A's spatial profile, and subtracted this spectrum from that of 
the nova.0 We found that Star A typically contributed 40% (65% on January 21) of the flux in 
the original extraction of Nova Vel 1993. In general (but, surprisingly, not on January 21), the 
subsequent analysis worked best on these fully decontaminated spectra of the nova, whose typical 
S/N ratio per final 0.75 A bin is 3.3 ± 0.4 in the continuum. 

3. RESULTS 

3.1. Unphased Average Spectra 

Two unphased average spectra of Nova Vel 1993 are shown in Figure 3, with that of 6 March 
1998 offset by 2.5 units. In the region of overlap, these resemble the spectra shown by Shahbaz et 
al. (1996). The strongest emission line is Ha, and there may be weak He I A5876 on 6 March, 
though the latter line partially coincides with the Na I D absorption. Despite being in quiescence, 
Nova Vel 1993 still exhibits variability; the equivalent width (EW) of the Ha emission line, which 
is unaffected by errors in flux calibration, was somewhat larger on March 6 (76 A) than on January 
21 (60 A). Note that this emission is considerably weaker than in GRO J0422+32 (EW « 250 A; 
Filippenko et al. 1995b), but comparable to that in GS 2000+25 (EW « 40 A; Filippenko et al. 
1995a) and Nova Oph 1977 (EW = 25-85 A; Filippenko et al. 1997). As in the previous objects 
we have studied, the Ha line has two peaks (Av ~ 1200 km s -1 ), more obvious in March than in 
January. 

3.2. Cross-Correlations 

Following the procedure discussed in Filippenko et al. (1995b), we employed the FXCOR 
package ("Release 9/13/93") in IRAF to cross-correlate the spectra of Nova Vel 1993 with the 
13 velocity standards (shown in Fig. 4, along with two others from our previous studies). The 
correlation was done over the ranges 5980-6270 A and 6320-6500 A to avoid the Ha and He I 



3 IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the Association of 
Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science Foundation. 

4 This extraction determined the flux scale of the subtracted spectrum, but its signal-to-noise (S/N) ratio was 
improved by using an extraction of the entire right half (rows 0-5) of Star A's spatial profile and scaling to the 
spectrum of rows 3-5. 
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A5876 emission lines, Na I D absorption, the 6270 A interstellar line, and poorly subtracted [O I] 
A6300 night-sky emission. In almost all cases a definitive correlation peak was obvious. Typical 
values of the Tonry & Davis (1979) significance threshold were quite high, R ;> 3, with a few as 
high as ~ 6. As in our previous studies, we adopted the FXCOR la uncertainties reduced by 
a factor of 2.77; the Fourier transform properties of Gaussians (the functions used in fitting the 
cross-correlation peaks) were used to determine that FXCOR overestimates the Tonry k, Davis 
uncertainties by this amount. 

The strongest formal correlation was obtained with BD+00 3090, which is officially listed as 
an MO V star (Upgren et al. 1972), but over our spectral range it also looks very similar to K7 V 
and K8 V stars; see Figure 4. Indeed, the correlation was insignificantly lower with K7 V and K8 V 
stars, but far inferior with K5 V and Ml V stars. (We did not observe K6 V and K9 V stars.) 
Thus, we conclude that the secondary star lies somewhere in the range K7 V through MO V, and 
possibly as early as K6 V. This is a slightly later spectral type than preferred by Shahbaz et al. 
(1996; late-G/early-K) and Delia Valle et al. (1997; later than G5-K0), though the former authors 
note that their derived mean density for the secondary (2.4 g cm" 3 ) suggests a K5 V star. Our 
phased-average spectrum of the secondary star (see below) supports the late-K classification. 

Using the radial velocities evaluated from the correlations with BD+00 3090 (corrected for 
the radial velocity of BD+00 3090 itself; 46.6 km s~\ Evans 1967), we conducted a non-linear 
least-squares fit (i.e., a \ 2 fit; Press et al. 1986, p. 521) to obtain the best cosine curve to match 
the data (Fig. 5). The four-parameter fit (zero point, semiamplitude, period, and phase) yielded a 
systemic velocity of 72 = 40.7 ± 4.0 km s" 1 , a semiamplitude of K2 = 475.4 ± 5.9 km s -1 , a period 
of P = 0.285206 + 0.0000014 d (6.84 hr), and a starting time (heliocentric Julian day) for the phase 
of T = HJD 2,450,835.0661 ± 0.0007, where T is defined as the point of maximum redshifted 
velocity. All the uncertainties are the formal la values derived from the \ 2 fit, but that of To may 
be an underestimate because the choice of the range in which to search for To is unconstrained by 
the data. A better measurement of the systemic velocity is 72 = 30.1 ± 5.0 km s _1 , the weighted 
average of values obtained with the nine best standard stars. Note that when Star A (Fig. 1) was 
cross-correlated with the velocity standards, its derived velocities were constant to within ±7 km 
s -1 (full range); thus, our results for Nova Vel 1993 are not an artifact of telescope position or 
instrument orientation, and our removal of the contamination by Star A (spectral type ~K1) was 
effective. 

It is interesting to examine the distribution of possible periods resulting from our series of 
observations. In Figure 6, which shows \ 2 versus trial period, there are several major groups of 
possible periods; the structure within each group results from ambiguities in the counting of cycles 
between widely separated epochs of observation (for example, from January through March 1998). 
Until we obtained the observations of 1999 January 21, the periods near 0.22 d and 0.285 d were 
equally probable, with those near 0.18 d and 0.4 d distinctly inferior. With the additional data, 
however, we are now reasonably confident that 0.285 d (6.84 hr) is correct. Based on observations 
of photometric modulation observed during the decline from outburst, Bailyn & Orosz (1995) 
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had speculated that the orbital period might be 1.6 hr or ~ 3 d, but these are clearly excluded 
by our radial velocity measurements. The possible presence of "superhumps" in the optical light 
curve obtained 4 months after the primary outburst led Delia Valle et al. (1997) to deduce an 
orbital period of 4 hr, closer yet still incorrect. Our period is essentially identical that found by 
Shahbaz et al. (1996; 6.86 ± 0.12 hr) from Gunn i?-band ellipsoidal modulations in quiescence. 
Note that the secondary in Nova Vel 1993 must be a dwarf; with an orbital period of only 6.84 hr, 
the compact primary would be inside a giant or subgiant. 

The formal reduced \ 2 {xt) f° r the best fit is 1.55 (17 points, 13 degrees of freedom). Given 
that the velocity uncertainties given by the Tonry & Davis (1979) method don't reflect external 
errors such as miscentering of the object in the slit, they could easily be too small. In our case, 
increasing them by only 25% yields a reduced xt °f 1-0- Moreover, had we assigned uncertainties 
to the adopted times (i.e., phases) of the observations, and performed the fit to simultaneously 
minimize residuals in both velocity and phase, we would have obtained a smaller value of xt- (The 
effective times of the observations can differ from the calculated midpoints due to variations in 
observing conditions — seeing, transparency, etc.) 

3.3. The Phased Average Spectrum 

Having determined the orbital parameters of the secondary star, we obtained its master 
"rest-frame" spectrum by averaging the 17 spectra after Doppler shifting each one to zero velocity. 
This represents a total integration time of ~5 hr, but of course the S/N ratio is lower than for a 
single 5-hr exposure (ignoring cosmic rays) due to the increase in readout noise. Note that the Ha 
and He I emission lines are smeared out by this process, since they are produced in the accretion 
disk around the compact primary star. Similarly, interstellar absorption lines become less distinct. 

Figure 7 shows the phased average spectrum of Nova Vel 1993, in comparison with its typical 
spectrum (obtained at 11:14 UT on 1998 Feb. 1), the unphased average around Ha, and the 
spectrum of the MO V velocity standard star. Scrutiny of the phased average spectrum reveals 
stellar absorption lines that are also present in the MO V star (e.g., the Ca + Fe blend at 6498 A), 
but those of Nova Vel 1993 are weak. While this could indicate that the spectral type of the 
secondary star is substantially earlier than MO, tests show that several other factors are likely to 
be more important: (1) rotational broadening, perhaps 50-100 km s -1 (e.g., Wade & Home 1988; 
Harlaftis et al. 1996); (2) orbital broadening, typically AK 2 T/P = 80-90 km s" 1 (Filippenko et 
al. 1995b); and, most importantly, (3) contamination by the featureless continuum of an accretion 
disk. 

We attempted to quantify the contribution of the accretion disk in Nova Vel 1993 by 
comparing its phased average spectrum with that of the M0 V star, broadening and diluting 
the latter by various amounts. A good match to the depths of the narrow absorption lines 
was found with accretion disk contamination fractions of 60%-70% of the total flux density at 
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~ 6300 A. Because of the uncertainties in the adopted parameters (broadening, spectral type) 
and the relatively low S/N ratio of the Nova Vel 1993 spectrum, we cannot confidently exclude 
contributions somewhat outside this range, but it is clear that the accretion disk dominates the 
spectrum even in quiescence. 

No absorption line of Li I A6708 is visible in our phased-average spectrum, to a 3a upper 
limit of ~ 0.1 A. An unexpectedly strong Li I line (EW = 0.25-0.48 A) is seen in the spectra of 
several X-ray novae (Martin et al. 1994, and references therein; Filippenko et al. 1995a), but not 
in others (such as Nova Oph 1977; Filippenko et al. 1997; Harlaftis et al. 1997). 

3.4. Ha Measurements 

To investigate the motion of the compact primary in Nova Vel 1993, we fit a Gaussian to 
the high- velocity wings (v ss 800-2000 km s _1 ) of the Ha line in each individual spectrum using 
the IRAF task SPECFIT. The regions of the fit (including the continuum) were 6400-6545 A 
and 6580-6710 A; the double-horned core of the line (Fig. 3) was excluded. The velocity of the 
Gaussian peak and its formal la uncertainty were adopted for each spectrum (Table 1). 

As shown in Figure 8, the derived velocity tends to be low in the first half of the orbit and 
high in the second half, perhaps suggesting periodic behavior. We used a least-squares fit to 
determine the Ha radial- velocity cosine curve, forcing the data to have the period found for the 
secondary star (0.285206 d) but allowing all other parameters to vary. The formal results (Fig. 8; 
xl = 1.34) are as follows: 71 = 4.6 ± 6.2 km s -1 , K\ = 65.3 ± 7.0 km s _1 , and a zero point in the 
phase of HJD 2, 450, 834.9686 ± 0.0093 d. The value of 71 is inconsistent with that of 72, but this 
is often the case in X-ray novae (e.g., Nova Oph 1977, Filippenko et al. 1997). 

Once again, the zero point in the phase is the maximum redshifted velocity; it implies that the 
compact object is 237° out of phase with the companion star, rather than the expected 180°. This 
is comparable to the offsets in GRO J0422+32 (253°; Filippenko et al. 1995b) and GS 2000+25 
(260°; Filippenko et al. 1995a), as well as in A0620-00 and Nova Mus 1991 (Orosz et al. 1994). 
To date, there is no satisfactory quantitative explanation for these distortions, but they suggest 
that the accretion disk often has a nonaxisymmetric distribution of surface brightness, noncircular 
velocities, or a warp. They also cast some doubt on the use of Ha radial-velocity curves to 
determine the motion and mass of the primary star. On the other hand, the mass ratios determined 
in this manner are frequently quite consistent with those obtained with independent techniques 
(e.g., A0620-00, Orosz et al. 1994; GS 2000+25, Harlaftis et al. 1996; Nova Oph 1977, Harlaftis 
et al. 1997; GRO J0422+32, Harlaftis et al. 1999). Thus, here we will cautiously adopt the ratio 
of semiamplitudes as an estimate of the mass ratio: q = M2/M1 = K1/K2 = 0.137 ± 0.015. 
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4. THE MASS OF THE COMPACT PRIMARY 

From the semiamplitude (K 2 = 475.4 ±5.9 km s" 1 ) and period (P = 0.285206 ± 0.00000138 d) 
of the radial velocity curve of the secondary star, we find a mass function /(Mi) = PK 2 /2ttG = 
3.17 ± 0.12 M . This corresponds to the absolute minimum mass of the compact primary, and 
it is close to the maximum gravitational mass of a slowly rotating neutron star (3.0-3.2 M ; see 
Chitre & Hartle 1976, and the discussion in Filippenko et al. 1995b). 

No evidence for eclipses is seen in the data of Bailyn Sz Orosz (1995), Shahbaz et al. (1996), or 
Delia Valle et al. (1997); similarly, we do not see significant variations in the apparent brightness 
of the secondary star over the orbital period. Hence, it is likely that the orbital inclination i <J 80°, 
and the relation /(Mi) = Mf sm 3 i/(M 1 + M 2 ) 2 then implies that M 1 ^ 4.2-4.4 M for nominal 
secondary-star masses of 0.5-0.65 M (M0-K6 V; Allen 1976). Even if the secondary is quite 
undermassive (e.g., M 2 = 0.3 M ), as in some X-ray binaries (e.g., van den Heuvel 1983), we 
derive Mi ^3.9 M . The primary star is therefore almost certainly a black hole rather than a 
neutron star. 

Adopting the mass ratio derived from the measured semiamplitude of the radial velocity 
curve of the primary (q = M 2 /M 1 = Ki/K 2 = 0.137 ± 0.015), we find Mi = 3.64-4.74 M if 
M 2 = 0.5-0.65 Mq. The constraints from q and the mass function yield Mi = 4.4 M and 
i ss 78° if we use a normal K7-K8 secondary (M 2 ss 0.6 Mq). Indeed, consistency cannot be 
achieved for M 2 <; 0.59 M if the maximum inclination estimate (80°) is correct. We conclude 
that the secondary star cannot be substantially undermassive, and that the orbital inclination is 
probably rather large, almost making Nova Vel 1993 an eclipsing binary. Our suggested inclination 
is significantly higher than the nominal value derived by Shahbaz et al. (1996) from their 
observed i?-band ellipsoidal modulations (i = 44° ± 7°). However, these authors admit that when 
contamination by light from the accretion disk is included, their allowed range for the inclination 
is 37°-80°. Of course, in view of the unexplained phase offset between the expected and observed 
Ha radial velocity curves (§ 3.4), it is also possible that our derived value of q is erroneous, thereby 
affecting our estimate of i. Further studies are needed to accurately determine the inclination. 

Recently, Bailyn et al. (1998) found that the distribution of masses of putative black holes in 
LMXBs is very strongly peaked at ~ 7 M Q , only V404 Cyg being a clear high-mass deviant. The 
number of objects in the sample is still quite small, but if our mass estimate for Nova Vel 1993 
is correct, then it appears to be a low-mass counterexample (Mi ~ 4.4 Mq). Another possible 
exception is GRO J0422+32 (Mi ~ 5 M ; Harlaftis et al. 1999). An independent measure of the 
mass ratio of Nova Vel 1993 from the rotational broadening of the secondary star's absorption 
lines (e.g., Nova Oph 1977, Harlaftis et al. 1997; GRO J0422+32, Harlaftis et al. 1999), together 
with better constraints on the inclination derived from near- infrared ellipsoidal modulations (e.g., 
GS 2000+25, Callanan et al. 1996), would provide a very useful check on our estimated mass for 
the primary star. 

When we calculate the effective Roche lobe radius (Rl) of the companion star from the 
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relation of Paczyriski (1971; see also Eggleton 1983) and Kepler's third law, we find that 
R L = 0.7 R Q if Mi ps 4.4 M Q and M 2 ps 0.6 M . This is only slightly larger than the expected 
radius of a typical K8 dwarf (R = 0.67 R Q ; Allen 1976). Thus, the secondary star may be just 
starting its evolution off the main sequence. 

5. OBSERVATIONS OF MXB 1659 29 IN QUIESCENCE 

As part of our effort to determine the mass functions of X-ray binaries, we observed the field 
of MXB 1659-29. This burst source was discovered in 1976 October with SAS 3 (Lewin et al. 
1976). It was initially considered unusual because of its very stable burst intervals (Lewin 1977) 
and apparent absence of constant emission, but such emission was found a year later (Lewin et al. 
1978; Share et al. 1978). An optical counterpart (now known as V2134 Oph) was discovered by 
Doxsey et al. (1979) at V = 18.3 mag; it was quite blue, and possibly exhibited emission lines of 
He II A4686 and C III/N III AA4640-4650. Figure 2 of Doxsey et al. (in which N is up and E to 
the left, although this isn't stated) shows [/-band and fi-band finder charts for the object, from 
images taken on 1978 May 30 and June 1 with the CTIO 4-m telescope. 

We obtained three dithered i?-band images (exposure times of 60, 60 and 30 s) of the 
MXB 1659-29 field on 1999 February 9 with LRIS/Keck-II, at airmass 1.8-1.9. These were 
bias-subtracted, flattened, registered, and combined in the usual manner. A small subset of the 
resulting crowded image (Galactic latitude 7.3°) is shown in Figure 9a; the FWHM of stars is 
measured to be ~ 0.9". Star A is close to the apparent position of V2134 Oph indicated in 
the relatively shallow charts published by Doxsey et al. (1979). However, seven LRIS/Keck-II 
long-slit spectra (PA = 160°) of this star, each with a typical exposure time of 900-1000 s, do 
not reveal any Ha emission characteristic of accretion disks. Moreover, cross-correlation of the 
individual spectra with the 13 velocity standards (G5 providing the best match) reveals no clear 
variability beyond the ±15 km s _1 level, and no systematic trend among consecutive exposures, 
casting further doubt on Star A as the secondary. Ha emission is also weak or absent in our noisy 
spectrum of Star B (Fig. 9a), and Star E has a spectral type of M. 

Shortly before the completion of this paper, MXB 1659-29 went into outburst again, after 
a hiatus of 21 years. During the interval 1999 April 2.06-3.47, the Wide Field Camera on 
BeppoSAX detected a transient X-ray source coincident with the position of MXB 1659-29 (in 't 
Zand et al. 1999). This object was confirmed on April 5.83-6.05 with RXTE (Markwardt et al. 
1999). Optical observations (Augusteijn, Freyhammer, & in 't Zand 1999) on April 3.41 revealed 
a bright new source (V = 18.3 ± 0.1) at that location, with a spectrum typical of LMXBs in an 
X-ray bright phase (emission lines of H I, He II, C III, and N III). An image (exposure time 30 s) 
obtained on April 18 with LRIS/Keck-II is shown in Figure 9b (stellar FWHM = 0.65"); the 
optical counterpart (R = 18.2 ± 0.05) is marked "F." This star is also visible at the center of the 
circle in Figure 9a, barely at the detection limit. 
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We measured the R magnitude of V2134 Oph in quiescence (Fig. 9a) with the technique of 
point-spread-function (PSF) fitting, where the PSF was determined iteratively by subtracting faint 
stars near the ones chosen for the PSF. The zero point was obtained from twilight-sky observations 
of PG1525-071A,C (Landolt 1992). Star F, which we identify with V2134 Oph in quiescence, has 
R = 23.6 ± 0.4 mag. If Star F is just a chance superposition, then the true optical counterpart 
is even fainter. Thus, any future attempts to obtain the mass function of MXB 1659-29 (V2134 
Oph) will be extremely difficult to perform! As an aid for future photometry of this object, we 
note that the final magnitudes for Stars A, B, C, D, and E are 19.6, 22.5, 23.1, 23.3, and 21.0, 
respectively. The la uncertainty is about 0.05 mag at the bright end (primarily due to the dearth 
of photometric standards) and perhaps 0.2 mag for stars at R 23. 

6. CONCLUSIONS 

Our observations of Nova Vel 1993 provide a definitive mass function of 3.17 ± 0.12 M , and 
a likely mass of around 4.4 M Q for the primary star. Thus, Nova Vel 1993 joins the small but 
growing list of secure Galactic black holes first identified as X-ray novae. However, its mass seems 
to be lower than that of other objects in its class, bridging the apparent gap between ~ 3 M 
(the theoretical maximum mass of a neutron star, though observed masses almost always yield 
M = 1.0—1.8 Mq; Thorsett et al. 1993) and ~ 7 M Q (the mass of most Galactic black holes in 
binary systems with well-determined parameters). It will be important to measure the mass ratio 
and orbital inclination of the system with techniques independent of the indirect ones used here, 
to confirm our estimates (q = 0.137 ± 0.015; i « 78°) and our derived mass. 

A similar study of MXB 1659-29 (V2134 Oph) eliminates several candidate stars as the 
optical counterpart. The recent new outburst, which occurred just prior to the submission of this 
paper, allows us to identify the quiescent nova at R = 23.6 ± 0.4 mag, unless this is an unrelated 
star superposed along the line of sight. 

Data presented herein were obtained at the W. M. Keck Observatory, which is operated as a 
scientific partnership among the California Institute of Technology, the University of California, 
and NASA. The Observatory was made possible by the generous financial support of the W. 
M. Keck Foundation. We thank T. Bida, R. Campbell, R. Goodrich, D. Lynn, G. Puniwai, H. 
Rodriguez, C. Sorensen, W. Wack, T. Williams, and other members of the Keck staff for their able 
assistance. We are also grateful to N. Vogt, A. C. Phillips, and D. C. Koo for obtaining a Keck 
image of MXB 1659-29 on April 18, after its recent outburst. This work was supported by the 
NSF through Grant AST-9417213 to A.V.F. 
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Table 1. Journal of Observations and Radial Velocities 



HJD a 


UT Date 


Exposure (s) 


Airmass b 


Phase c 


v 2 (km s- 1 ) 6 


vi (km s^ 1 ) 6 


838.96576 


1998 Jan 25 


900 


2.53 


0.67402 


-236.7 ±21.7 


123.0 ±62.3 


845.95655 


1998 Feb 01 


1200 


2.46 


0.18531 


219.4 ± 18.9 


-36.10 ±21.7 


845.97085 


1998 Feb 01 


1200 


2.39 


0.23545 


52.09 ± 17.1 


-21.29 ± 25.1 


877.85456 


1998 Mar 05 


900 


2.58 


0.026917 


448.5 ± 13.4 


-61.26 ±28.3 


877.86708 


1998 Mar 05 


1200 


2.48 


0.070815 


430.6 ± 15.0 


-79.79 ± 19.7 


877.88132 


1998 Mar 05 


1200 


2.41 


0.12074 


328.9 ± 8.62 


-34.98 ± 18.8 


877.89554 


1998 Mar 05 


1200 


2.36 


0.17060 


248.4 ± 17.5 


-51.49 ± 17.6 


877.90868 


1998 Mar 05 


1000 


2.35 


0.21667 


123.3 ± 15.3 


-71.38 ±21.1 


877.92060 


1998 Mar 05 


1000 


2.35 


0.25847 


-61.19 ±15.5 


-86.64 ±21.0 


878.85228 


1998 Mar 06 


1100 


2.58 


0.52515 


-476.1 ± 12.0 


75.89 ±25.1 


878.89168 


1998 Mar 06 


1100 


2.37 


0.66329 


-229.7 ± 17.2 


63.71 ±22.3 


878.90475 


1998 Mar 06 


1100 


2.35 


0.70912 


-120.3 ± 15.8 


40.81 ± 20.7 


878.91785 


1998 Mar 06 


1100 


2.36 


0.75505 


11.61 ±21.8 


86.20 ±18.5 


1199.97900 


1999 Jan 21 


1200 


2.52 


0.46778 


-476.9 ± 12.0 


14.13 ±41.8 


1200.00574 


1999 Jan 21 


1100 


2.38 


0.56154 


-451.1 ±8.61 


27.51 ±50.5 


1200.01881 


1999 Jan 21 


1100 


2.35 


0.60736 


-370.1 ± 12.8 


31.60 ±66.0 


1200.03255 


1999 Jan 21 


1200 


2.35 


0.65554 


-272.4 ±27.6 


14.54 ± 63.6 



a HJD-2,450,000 at midpoint of exposure. 
b Airmass at midpoint of exposure. 
c Using P = 0.285206 d and T = 2,450,835.0661. 
d Secondary star radial velocity. 

e Ha centroid radial velocity, from fit to emission-line wings. 
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Fig. 1. — i?-band image of Nova Vel 1993, obtained on 1998 January 25. The now-quiescent nova 
is about 1.7" NNW of Star A. 

Fig. 2. — Spatial cut of counts along the slit in a typical two-dimensional CCD spectrum of Nova 
Vel 1993, computed from the average of 10 columns centered on the wavelength of Ha. The nova 
is in the wing of Star A. Other stars correspond to those shown in Fig. 1. 

Fig. 3.— Typical spectra of Nova Vel 1993, obtained on 1998 March 6 and 1999 January 21. Note 
the change in the shape and strength of the Ha emission line. The top spectrum is offset by +2.5 
units for clarity. 

Fig. 4. — Spectra of 15 radial velocity standards, arranged from early to late types. 

Fig. 5. — Radial-velocity curve of the secondary star in Nova Vel 1993, derived from cross- 
correlations of individual spectra with those of the MO V star BD+00 3090. Two cycles are shown 
for clarity. The radial velocity of the comparison star, 46.6 km s _1 , has been added to all values. 
The orbital period is 0.285206 ± 0.00000138 d. Formal velocity error bars are la. 

Fig. 6. — Value of x 2 vs. trial period for Nova Vel 1993, using the 17 spectra discussed in this 
paper. Only the lowest values of x 2 are shown; x 2 ~ 10, 000 at typical periods near 0.35 d. 

Fig. 7. — Spectra of Nova Vel 1993 compared with a scaled spectrum of a M0 V velocity standard. 
Constants have been added to the top three spectra for clarity. The "phased average" spectrum 
in the rest frame of the companion star was obtained by combining the 17 spectra after Doppler 
shifting each one to zero velocity. The narrow Ha absorption in the unphased average spectrum 
(top) must be produced by gas associated with the accretion disk, not the secondary star; otherwise, 
it would appear smeared out. 

Fig. 8. — Radial-velocity curve of the centroid of the Ha emission line, determined by fitting a 
Gaussian to the high-velocity wings. The period was forced to be that of the secondary star, but 
all other parameters were allowed to vary. Formal error bars are la. 

Fig. 9.— (a) iZ-band image of the field of MXB 1659-29 (V2134 Oph), obtained on 1999 February 
9. Various stars discussed in the text are marked, (b) i?-band image obtained on 1999 April 18. 
The nova in outburst is positionally coincident with Star F in (a). 
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